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Abstract. We study how the volatility, node- or link-based, affects the evolution of social networks in 
simple models. The model describes the competition between order - promoted by the efforts of agents 
to coordinate - and disorder induced by volatility in the underlying social network. We find that when 
volatility affects mostly the decay of links, the model exhibit a sharp transition between an ordered phase 
with a dense network and a disordered phase with a sparse network. When volatility is mostly node-based, 
instead, only the symmetric (disordered) phase exists These two regimes are separated by a second order 
phase transition of unusual type, characterized by an order parameter critical exponent /3 = + . We argue 
that node volatility has the same effect in a broader class of models, and provide numerical evidence in 
this direction. 



PACS. 89.65. 
Networks 



1 Introduction 



s Social and economic systems - 05.70.Fh Phase transitions: general studies - 64.60.aq 



The competition between order and disorder is by no means 
restricted to physics. Also economies and societies - as 
systems of many interacting individuals - organize them- 
selves in different (macroscopic) states, with different de- 
grees of order - informally interpreted as coordination on 
social norms, compliance with laws or conventions pQ. Be- 
sides all its inherent complexity, one important element of 
additional richness is that the relation between the degree 
of order in a society and the cohesion of the underlying 
social network is not unidirectional as in physics, where 
the topology of interactions is fixed [3] . Rather the degree 
of order in a society influences in important ways the den- 
sity and topology of interactions. The interplay between 
network's dynamics and collective behavior is important 
in many phenomena, ranging from informal contacts in 
labour market [2] and peer effects in promoting (anti- 
social behaviors [6; to inter-firm agreement for R&D [5]. 
The structure of the networks involved in these phenom- 
ena is dynamically shaped by incentives of agents (nodes) , 
be they individuals or organizations, who establish bilat- 
eral interactions (links) when profitable. 

In addition, this interplay typically takes place in a 
volatile environment. That is, the favourable circumstances 
that led at same point to the formation of a particular link 
may later on deteriorate, causing the removal or rewiring 
of that link. This combination of factors raises a number 
of interesting issues in statistical physics, as the collec- 
tive behavior - of e.g. processes of ordering [13], opinion 
spreading, [4] and reaction diffusion [5] - may radically 



change when they are coupled to the dynamics of the net- 
work they are defined on. 

We shall here focus on the stylized mathematical de- 
scription of this generic phenomena given in Ref. [13j : Here 
the feedback between nodes and networks dynamics, arises 
from assuming that the formation of a link requires some 
sort of similarity or proximity of the two parties. This cap- 
tures different situation: For example, in cases where trust 
is essential in the establishment of new relationships (e.g. 
in crime or trade networks), linking may be facilitated 
by common acquaintances or by the existence of a chain 
of acquaintances joining the two parties. In other cases 
(e.g. in R&D or scientific networks) a common language, 
methodology, or comparable level of technical competence 
may be required for the link to be feasible or fruitful to 
both parties. 

This class of models reveals a generic behavior char- 
acterized by a discontinuous transition from an uncoor- 
dinated state characterized by a sparse network, to a co- 
ordinated state on a dense network. As discussed in Ref. 
[13] , this agrees with anecdotical evidence which can be 
summarized as follows: 

(i) Sharp transitions. Observation on the spread of so- 
cial pathologies , the growth of research collaborations, 
both scientific [7] and industrial [8 ], suggest that networks 
can shift from a sparse to connected state in short time 
spans. 

(ii) Resilience. Once a transition to a highly connected 
network has taken place, this setup can survive even to a 
reversion to unfavorable conditions, e.g. the thriving per- 
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formance of Silicon Valley during the computer industy 
crisis of the 1980sj9], or the recent development of open- 
sourcesoftware, sustained against large odds, thanks to a 
dense web of collaboration and trust [ID]. 

(iii) Equilibrium coexistence. Under apparently similar 
environmental conditions, social networks can be found 
both in a dense or sparse state. A good illustration is pro- 
vided, e.g. by the dual experience of poor neighborhoods in 
large cities, where neither poverty nor other socioeconomic 
conditions alone can explain wheter there is a degradation 
in a ghetto with rampant social problems [6] . 

This paper focuses on analyzing the effect of node 
volatility in these simple models. Indeed, the effect of 
volatility is limited to link removal in Ref. 13J, but the 
turnover of agents (i.e. node removal and arrival) may 
be an important factor in many real systems. Our focus 
here will be mostly on the statistical phenomenon, than 
on its interpretation in socio-economic terms. Indeed we 
find that the introduction of node volatility brings in a 
qualitative change, which can be described as a continu- 
ous phase transition with unusual critical properties. In 
order to show this, we concentrate on the simplest model 
for which a full analytic treatment is possible. In the con- 
cluding section, we argue that this qualitative change is 
expected in a wider class of model, and it can have much 
stronger effects. 



2 model 

Our model is a variant of one in the general class of [13]. 
It reproduces in a stylized manner the mechanisms of co- 
evolution in social networks aforementioned, and shows 
some common elements of the observed phenomenology 
(i.e. sharp transitions, resilience, equilibrium coexistence); 
in particular, it also shows how the alternative assump- 
tions of link or node based volatility have profound effects 
on the dynamics of network formation. The model de- 
scribes N agents sitting on the nodes of a network, each 
of which is characterized by a variable o~i which represents 
the social norm (convention or technological standard) 
adopted by agent i. There are q possible social norms, i.e. 
&i G {1, . . . , q}. In terms of statistical physics, the model 
can be thought of as a q state Potts model defined on a 
graph of N nodes, that evolves in a coupled fashion to the 
dynamics of the system. The rules of the dynamics are the 
followings: 

— each node attempts to establish a new link with a ran- 
domly chosen node at rate 77/2. The link is established 
only if the two nodes have the same color 

— links are destroyed at rate 1 

— all the links of a node are destroyed with rate a 

— the color of a node is updated with rate v to the color 
of any of its neighbors, unless the node is isolated. In 
the latter case the nodes takes a random color. 

So, the parameter a interpolates between two kinds 
of volatility. For a = volatility only affects links and 
for a > 1 it mostly affects nodes. With respect to the 



parameters introduced in Ref. [13] . we observe that the 
link decay rate has been set to A = 1 and that, for no- 
tation convenience, we scaled by a factor 2 the link cre- 
ation rate r\. As observed in Ref. [13], the color update 
rule is effective only for isolated nodes, in the long run, 
and in that case the color is drawn at random. Since only 
links between same type agents are created, after a tran- 
sient all nodes are either isolated, or connected to nodes 
of the same color. Therefore the particular way in which 
the neighbor is chosen is immaterial. For example, both 
a majority rule (most frequent color) or a voter-type rule 
(random neighbor) would give the same dynamics. The 
model can be generalized to a probabilistic update rule 
for the colors introducing a finite temperature T (see [13]). 
Results do not change considerably as long as T is small 
enough, so we shall confine ourselves to the T = case. 




Fig. 1. Our model: links are formed with rate n between nodes 
having the same color, links are destroyed with rate 1, all the 
links of a node desappear with rate a and colors of isolated 
nodes are randomly updated with rate v. 



Ref. [13] has shown that for a = 0, the system shows an 
hysteretic transition in 7/ from a symmetric to an asym- 
metric state. The symmetric state is characterized by a 
sparse network, with average degree (k) < 1, with a sym- 
metric distribution of colors. In the asymmetric state, in- 
stead, a dense network with (k) > 1 arises, along with 
a dominant color, which is adopted by agents more fre- 
quently than the others. In this sense, the model shows 
how order and disorder are intimately related with the 
dynamics of the social network in a volatile environment. 

In what follows, we solve the model in the stationary 
state for N — > 00, for all the values of a. We find that the 
a = behavior is generic for all a < 1, but the transi- 
tion is softened as a increases. For a > 1 instead we show 
that the system is always in the symmetric phase. Hence, 
in terms of statistical mechanics, a = 1 is a second or- 
der critical point separating a phase with spontaneously 
broken symmetry from a symmetric phase. 
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3 Theory 

If wc call rifc ]0 . the density of nodes with k links and color 
a = 1, . . . , q we have the following rate equations: 

flk,a = (k + l)nk+l,a - kl%k,a - a"-fe, CT + 

+ X a (n,k-l,a -nk,a) (1) 

«o,o- = a(n a - n .<j) + ni jCT - x a n k: a + 



- ^2 (no t< T' ~ n . a ) 



Q 



(2) 



where, for future convenience, we have introduced the dy- 
namical variables 






(3) 



aking the sum over all k of these equations and multiplying 
by ?7 we find 



— > _, {n .cr> - n 0t a) 
q z -^ 



(4) 



which implies that, in the stationary state, each compo- 
nent has the same fraction 77,0,0- = n$/q of disconnected 
(k = 0) nodes. This is equivalent to a detailed balance 
condition for the density of the different components. 

It is straightforward to derive an equations for the 
characteristic functions Tr a (s) of the degree distribution 
p a (k) — nk l(7 /J2q n <i-<? °f the component a. In the sta- 
tionary state this reads: 

(1 - s)— — = [a + x a (l - s)]TT a (s) - a. (5) 

as 

The stationary solution is found by direct integration: 

7r CT (s) = a J dzz a - 1 e- x ' {1 - s){1 - z) (6) 

Jo 

It is easy to see that this interpolates between a Pois- 
son distribution, 7r ff (s) = e x <?( s - 1 )/ x for a — > 0, which 
coincides with the result of Ref. [13] . and an exponen- 
tial distribution 7r CT (s) = A/ [A + x a (s — 1)] for a — » 00. 
Notice also that the average degree in component a is 
(k) a = TTo-(l) = x a /(l + a). This is precisely what one 
expects from balance of link creation and destruction of 
links in component a. 

Observing that ^(0) = ij^ 22 - = ^^ we find an equa- 
tion for x a in the stationary state, which reads 



G a {x a )=ax a / duu a - 1 e x ' {u - 1) 







rjno 



(7) 



Notice that the r.h.s. of Eq. (J7]) is independent of a. The 
variables x a are determined by Eq. ([7]) and the normal- 
ization condition, which takes the form 



1 

(7=1 



(8) 



The properties of the solutions of Eqs. (|7|8p depend on 
the behavior of the function G a (x), which are discussed 
in the appendix, and can be classified in symmetric and 
asymmetric solutions. 



3.1 a > 1: The symmetric solution 

For a > 1 the function G a (x) is a monotone increasing 
function (see appendix), Hence Eq. ([7]) has a single solu- 
tion and Eq. ([8]) implies that x a = i]/q for all components 
a. Hence Eq. (UJ yields the total fraction of disconnected 
nodes 

n = -G a (r)/q) 
V 

as a function of the parameters q, a and n. We can ana- 
lyze the stability of the symmetric solution recalling that 
n o,cr = "nG a (x a ). Then Eq. ^ becomes a dynamical equa- 
tion for x„ 



1 



9 






\G a (x a i) — G a \x a )\ 



(9) 



Linear stability of the symmetric solution is addressed by 
setting x a — rj/q + e a , with X)o- e °" = 0- Then to linear 
order 



11 
-( 

(1 



' = -G'^/q) Y, t £ -' - e -l = -vG'Jn/fr*- ( 10 ) 



Hence, as long as G a (x) is an increasing function of x, the 
symmetric solution is stable. This is always the case for 
a > 1, as we shall see, it fails to hold for a < 1. 



3.2 a < 1: The asymmetric solution 

For a < 1 the symmetric solution still exists. However the 
function G a (x) now has a maximum for some xq{cx) (see 
the appendix) and G a (x) — > a from above as 1 — ► 00. 
Therefore the symmetric solution becomes unstable when 
rj > r/ + where 

r/ + = ?a;o(a) (11) 

because beyond that point G' a (rj/q) < 0. 

The occurrence of a maximum in G a also implies that 
Eq. ((7j admits solutions with x a — X- < xq((x) for some 
cr's and x a = x + > x§{ol) for the other components. Since 
x a is related to the density of a component a, we shall call 
a component dense if ie^ = x + and diluted if x a = X-. 
As in Ref. [13], all solutions with more than one dense 
component are unstable. Indeed, by the same argument 
used to analyze the stability of the symmetric solution, a 
perturbation with e CT = for all diluted components would 
grow as e CT = — vG' a {x + )e a on all dense components. These 
unstable modes correspond to density fluctuations across 
dense components and are clearly absent in the solution 
with one only dense component. These are the asymmetric 
solutions we shall focus on in what follows. There are q 
of them, depending on which color is associated with the 
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dense component. The variables x± are determined by the 
system of equations 

G a (x+) = G a {x-) (12) 

x + + (q— l)x- = r) (13) 

This solution is shown in Fig. [21 Actually, of the two asym- 




Fig. 2. Solutions x± as a function of r\ for q — 10 and a — 
0.4. The dashed line x = rj/q separating the two curves is the 
symmetric solution. 



metric solutions the one with x+ decreasing with 77 is 
clearly unphysical as this would have a connected com- 
ponent with an average degree {k) a = x + /{l + a) which 
decreases with the rate 77 with which links are formed. As 
in ref. [13 , it is easy to see that only solutions with x+ 
increasing in 77 are stable. Indeed, regarding 77 and x_ as 
functions of x+ in Eq. (|13|) we find 



di] 
dxj. 



1 + (q - 1 



t dx- _ G' a (x-) + (q-l)G' a (x + ) 

dx + G' a (x-) 



Consider perturbations of the form x a = x + + e for the 
dense component and x a = X- — e/(q — 1) for the others. 
Then by a derivation analogous to that leading to Eq. 
(fTOl). we find 



e = - [G' a (x.) + (q- l)G' a (x + )} e = - V -G' a {x_)-p-e. 

q q dx+ 

Given that G' a {x-) > 0, this implies that on solutions 
with x+ decreasing with 77, the perturbation e grows un- 
boundedly. 

The asymmetric solution ceases to exist for rj < i]-\j. 
In the region 77 G [v-j 7 1+] both the symmetric and the 
asymmetric solutions co-exist. The coexistence region, in 
the a, r\ plane is reported in Fig. 2J 

The practical relevance of the results derived so far is 
best discussed introducing an order parameter 



(14) 



V 



which is the difference in the density of the dense and 
diluted components. This vanishes in the symmetric phase 



We note, in passing, that the condition dri/dx^ 
vides an equation which allows to determine r\- . 



pro- 



and is non-zero in the asymmetric one. In Fig. [3] where we 
report the behavior of the average degree of the network 



(k) = y^rtfc.crfc 



i] 1 + (q — l)m 2 
+ a q 



(15) 



Fig. shows that as 77 sweeps through the coexistence 
region the system undergoes an hysteresis loop: the de- 
gree jumps from low to high values at 77+ as 77 is increased 
whereas when 77 decreases from large values, the network 
collapses back to the symmetric phase when r\- is crossed. 
In the case a = [13j . the symmetric phase is character- 
ized by sparse networks, with a vanishing giant compo- 
nent. This is no more true when < a < 1, specially close 
to ?7 + q Numerical simulations fully confirm this picture, 
even though for finite systems the symmetric (asymmet- 
ric) phase is meta-stable close to 77+ (7/-) and therefore 
the transition occurs for lower (larger) values of 77. 




Fig. 3. Mean degree < k > as a function of 77/A for a system 
with q = 10 colors, for a — and 0.2, simulations are for 
systems of 1000 nodes. 




Fig. 4. Phase diagram for q = 10. The symmetric phase ex- 
tends below and to the right of the (full) line rj- (a) whereas 
above the (dashed) line 77+ (a) only the asymmetric phase is 
stable. The coexistence region, where both phases are stable, 
is delimited by the two curves. 



2 Indeed the condition for the presence of a giant component 
is {k(k — l)) a > (k)cr which, by a straightforward calculation, 
reads 77 > g(l + a/2). At the critical point 77+ = qxo(a) this 
reads xo(a) > 1 + ct/2 which holds true for all a > 0. 
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3.3 The critical region: awl 

The behavior of the order parameter m on the critical 
lines which confine the coexistence region is shown in Fig. 
[5j This shows that the transition is continuous but with 
a peculiar critical behavior. In order to shed light on this, 
the appendix shows that, asymptotically for a ~ 1 

■q + (a) =qx {a) ~ -glog(l - a) + clog | log(l - a)| (16) 

with c > a constant. A detailed asymptotic analysis of 
the limit a — > 1 (see appendix) also shows that 



m ~ l/x (a) ~ | log(l - a) | 



(17) 



In terms of the usual description of critical phenomena, 
this model is consistent with an 



1 



\P 



where m 

exponent j3 = + . Indeed, the singular behavior of m is 

very close to that of a first order phase transition. 




way: while the network promotes similarity or proximity 
between nodes, proximity or similarity enhances link for- 
mation. This feedback allows the system to cope with envi- 
ronmental volatility, which acts removing links at a con- 
stant rate. Interestingly, the emergence of an "ordered" 
state plays a key role in this evolutionary struggle. 

We believe the general findings discussed here will ex- 
tend to the general class of models of Ref . [T3] . In partic- 
ular, we expect the phase transition to be blurred by the 
effect of node volatility and to disappear when the latter 
exceed a particular threshold. 

Actually, Fig. [6] shows that this is the case even for 
the model of Ref. [T2] . This is a model where link creation 
occurs either by long distance search at rate rj (as in the 
model discussed here) or through local search (on second 
neighbors) at rate £ . Again links decay at unit rate. We 
refer the interested reader to Ref. [H] for further details, 
for the present discussion let it suffice to say that the ef- 
fects of (link) volatility are contrasted by the creation of 
a dense network with small-world features (a somewhat 
similar model with node volatility has been considered in 
Ref. [H]). Fig. [5] shows that the effects of node volatility 
are very strong. Indeed, even a very small a reduces con- 
siderably the size of the coexistence region and the value 
a c at which the latter disappears is also relatively small. 



<k> T| = 0.01 



Fig. 5. Order parameter m on the boundary of the coexistence 
region ry_ (full line) and r/+ (dashed line) as a function of a for 
q = 10. 



a = 0.01 0.1 



4 Conclusions 



a a 8 a ' , , H e ' g a 



The introduction of node volatility, in the simple model 
discussed here, makes the transition from a symmetric 
(disordered) diluted network to an asymmetric (ordered) 
dense network less sharp. Indeed when node volatility dom- 
inates, i.e. when the number of links lost per unit time by 
node decay outnumber those lost from link decay (a > 1), 
the transition disappears altogether, and the symmetric 
(disordered) state prevails. The phenomenology is strongly 
reminiscent of that of first order phase transitions (e.g. 
liquid-gas or paramagnet-ferromagnet) though the criti- 
cal behavior is highly non-trivial. 

The virtue of the particular model studied is that it al- 
lows a detailed analytic approach which allows one to gain 
insight on all aspects of its behavior. This model belongs 
to a general class of models which embody a generic feed- 
back mechanism between the nodes and the network they 
are embedded in, which can be expressed in the following 



<k> 




Fig. 6. Mean degree as a function of the rate £ of formation 
of links with neighbours of neighbours, for N — 1000 (A = 1). 
Right: r\ = 0.01, Left: r) = 0.1. 
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These results suggest that node volatility is indeed a 
relevant effect in the co-evolution of socio-economic net- 
works, as it may affect in dramatic ways the ability of the 
system to reach a dense and/or coordinated state. 



Appendix 

The function G a {x) can be written as 

G a {x) = a I dz I 1 J 

For a > 1 we have 

dG a a(a — 1 



dx x 1 

Notice also that 

d G a (x) 
dx x 



dzz ( 1 



>0 



-a / du« (1 - w) a_1 e~ ux 



i.e. no = G a (rj / 'q) / '(rj / 'q) in the symmetric solution is a 
decreasing function of r\jq. In addition G a (x) ~ x for 
i<1, i.e. no — ► 1. 

For a = 1 — e we can approximate 



1 - e log 1 



G a (x) ~ (1-e) / du 
Jo 



= (l-e)(l-e-* + e(E i (x)- 1 )) 

where Ei(x) is the exponential integral function, and for 
e — > we have x± — > oo and Ei(x) ~ — so 

G a ~(l-e)(l-e- x + e/a;) 

We have at the critical point .t+ + (<j — l)x_ = gxo where 
rro is such that G' a (x^) — so e^-j- = 1 and xo ~ — log e + 

2 log | loge| by definition m = x °~ ;r ~ and x_ = xo(l — m) 
x + = xo(l + (g— l)m) so from G a (x + ) = G Q (x_) we have 

e~ x+ - e/x+ ~ e _:E " - e/a;_ 

and then 

qmXQ _mi„/, 5 -ijmio\ 



(l-m)(l + (g-l)m) 

we have mx — * c where c is given by 
c = e c (l - e- ?c )/g 



(1-e- 
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